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Ab  Initio  Molecular  Dynamics  Simulation  of  the  Amorphous 
Structure  of  Ca-Mg-Cu  and  Ca-Mg-Zn  Alloys 


O.N.  SENKOV  and  Y.Q.  CHENG 

The  atomic  and  electronic  structures  of  several  Ca-Mg-TM  amorphous  alloys  (TM  is  Cu  or  Zn) 
have  been  analyzed  using  ab  initio  molecular  dynamics  simulation  and  neutron  diffraction. 
Partial  pair  distribution  functions  have  been  produced  and  the  pair  bond  distances  and  partial 
coordination  numbers  have  been  reported  for  these  alloys.  Similarities  and  differences  in  the 
amorphous  structures  of  the  Ca-Mg-Cu  and  Ca-Mg-Zn  alloys  have  been  discussed.  Strong 
interactions  between  Ca-Cu,  Mg-Cu  and  Ca-Zn  atom  pairs  rooted  from  the  orbital  hybrid¬ 
ization  of  the  s-p-d  electrons  have  been  recognized  to  result  in  noticeable  shortening  of 
respective  atom  pair  bond  distances  and  pronounced  chemical  short  range  ordering  near  the  TM 
atoms.  Voronoi  tessellation  analysis  has  shown  that  the  polytetrahedral-type  clusters  and  five- 
coordinated  atom  pairs  dominate  in  the  amorphous  structures,  which  indicates  that  tetrahedra 
and  pentagonal  bi-pyramids  are  the  main  building  blocks  in  these  amorphous  alloys. 
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I.  INTRODUCTION 

Ca-Mg-BASED  bulk  metallic  glasses  (BMGs) 
have  unique  properties.  They  are  built  on  two  simple 
metals,  Ca  and  Mg,  while  all  other  BMGs  are  transition- 
metal  based  alloys.11’21  Many  Ca-Mg-based  alloys  have 
very  good  glass  forming  ability  (GFA),  and  even  their 
ternary  alloys,  e.g.,  Ca-Mg-Zn  or  Ca-Mg-Cu,  can  be 
made  fully  amorphous  when  they  are  cast  in  -8  to 
10  mm  thick  plates[3  81  or  up  to  15  mm  diameter 
rods.19'101  These  glasses  have  very  low  Young’s  and 
shear  moduli  that  are  comparable  with  the  moduli  of 
human  bones/11  131  low  density  (1.6  to  2.8  g/cm3)/11  and 
strong  relaxation  dynamics  of  the  super-cooled  li¬ 
quid/141  GFA,  as  well  as  physical  and  mechanical 
properties  of  metallic  glasses,  is  believed  to  depend  on 
the  type  of  short  range  order  (SRO)  and  medium  range 
order  (MRO)  of  alloying  elements  in  the  amorphous 
structure.  It  is  therefore  interesting  to  analyze  the 
amorphous  structure  of  Ca-Mg-based  BMGs. 

Structural  analysis  of  metal-metal  BMGs  is  mainly 
focused  on  transition  metal  glasses  such  as  Zr-Cu/15  181 
Zr-Pt/191  and  Zr-Cu-Al/20,2'1  The  results  indicate  that  the 
packing  of  atoms  in  these  materials  is  not  random,  but  is 
strongly  influenced  by  chemical  interactions.  In  particu¬ 
lar,  icosahedral  SRO  has  been  identified  and  correlated  to 
their  good  GFA/221  On  the  other  hand,  the  atomic 
structure  of  Ca-Mg-Zn  BMGs,  identified  with  the  use  of 
X-ray  and  neutron  diffraction  and  Reverse  Monte  Carlo 
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(RMC)  simulation,  shows  no  icosahedral  SRO/2’1 
Instead,  five-fold  bonds  in  the  form  of  pentagonal 
bi-pyramids  have  been  found  to  be  the  most  populous 
structural  units  in  these  BMGs/2’1  A  large  fraction  of  five¬ 
fold  bonds  and  the  lack  of  icosahedral  SRO  has  also  been 
found  in  the  molecular  dynamic  (MD)  simulated  amor¬ 
phous  structure  of  Mg-Cu  alloys/241 

In  the  present  work,  the  atomic  structures  of  several 
Ca-Mg-Cu  and  Ca-Mg-Zn  BMGs  were  simulated  with 
the  use  of  ab  initio  (quantum)  molecular  dynamics 
(QMD)  simulation.  The  simulated  structures  were  then 
used  to  calculate  partial  (PRDF)  and  total  (RDF)  radial 
distribution  functions,  pair  bond  distances  and  partial 
and  total  coordination  numbers.  The  QMD-simulated 
results  were  validated  through  the  experimentally  deter¬ 
mined  RDFs. 


II.  EXPERIMENTAL  PROCEDURES 

QMD  simulation  of  the  amorphous  structures  of  the 
Ca60Mg15Zn25,  Ca6oMg25Zn15,  Ca6oMgi5Cu25,  Ca60Mg25 
Cu15,  Ca5oMg25Cu25,  and  Ca4oMg25Cu35  alloys  was  con¬ 
ducted  using  the  Vienna  ab  initio  Simulation  Package 
(VASP)/25]  The  Projector  Augmented -Wave  (PAW) 
method126’271  and  Perdew-Burke-Ernzerhof  (PBE)  func- 
tional128,291  were  used  to  describe  interacting  valence 
electrons.  Cubic  simulation  boxes  with  periodic  boundary 
conditions  contained  200  atoms,  and  the  box  size  was 
determined  by  the  density  of  the  modeled  material.  The 
simulation  was  conducted  at  F  point  only,  with  a  first  order 
Methfessel-Paxton  smearing  function  (a  =  0.2  eV).  The 
liquid  was  equilibrated  at  1000  K  (727  °C),  then  was 
quenched  to  300  K  (27  °C)  in  100  K  (100  °C)  temperature 
steps,  equilibrating  at  each  temperature  for  3000  steps. 
Another  3000  consecutive  configurations  were  further 
sampled  at  300  K  (27  °C)  with  a  time  step  of  5  fs,  and 
averaged  PRDFs,  g//MD(r),  were  calculated  from  these 
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Table  1.  Density  (in  g/cm3  and  Atoms/ A3)  of  Ca-Mg-TM 
(TM  is  Cu  or  Zn)  Amorphous  Alloys  Studied  in  This  Work 


Alloy 

Density 

g/cm3 

Atoms/A3 

Ca4oMg25CU35 

2.936  ±  0.007 

0.03987 

Ca5oMg25Cu25 

2.439  ±  0.003 

0.03497 

Ca6oMg25Cui5 

2.039  ±  0.003 

0.03097 

Ca60Mg15Cu25 

2.367  ±  0.004 

0.03271 

Ca6oMg25Zn15 

2.004  ±  0.004 

0.03023 

Ca60Mg15Zn25 

2.289  ±  0.007 

0.03130 

configurations  using  Eq.  [1],  thus  taking  into  account 
thermal  vibrations  and  assuming  no  structural  relaxation 
within  15  ps. 


gy(r) 


1 


driijir) 


4nr2Cjp0  dr 


to  reduce  termination  ripples  in  the  Fourier  transform 
(Eq.  [3]). 

Using  Eq.  [1],  the  partial  coordination  number  of  an 
element  j  in  the  first  coordination  shell  of  the  element  i 
was  determined: 


*  max 

CN  =  4nCjp0  J  gif2 dr  [4] 

^min 


Here  rmin  and  rmax  are  the  positions  of  the  start  and 
end  of  the  first  peak  in  the  respective  gij(f). 

The  simulated  structures  were  statistically  analyzed 
using  pair  and  three-body  correlation  functions,  Voro- 
noi  tessellation  and  nearest  neighbor  approaches,  which 
allowed  calculation  of  local  structural  features  (SRO 
and  MRO)  such  as  total  and  partial  CN,  type  and 
distribution  of  characteristic  coordination  polyhedra, 
atomic  volume,  and  local  packing  fraction. 


Here  day  are  the  number  of  elements  of  type  j  between 
distances  r  and  r  +  dr  from  an  element  of  type  i,  Cj  is  the 
atomic  fraction  of  the  element  of  type  j,  and  pQ  is  the  alloy 
density.  The  QMD  models  provided  a  complete  set  of  six 
simulated  PRDFs  for  each  of  the  studied  alloys,  which 
well  agreed  with  the  experimental  diffraction  data.  Due  to 
the  limited  box  size  and  the  periodic  boundary  condi¬ 
tions,  g^MD(r)  were  calculated  at  r  <  9  A. 

Fully  amorphous  samples  of  the  studied  Ca-Mg-Cu 
and  Ca-Mg-Zn  alloys  were  prepared  by  melt-spinning  in 
the  form  of  ribbons  as  described  in  Reference  23.  The 
densities,  p0,  of  the  amorphous  alloys  were  measured 
with  a  helium  pycnometer  AccuPyc  1330  VI. 03  and  the 
values  (in  g/cm3  and  atoms/A3)  are  given  in  Table  I. 
Neutron  diffraction  experiments  were  conducted  at 
room  temperature  under  vacuum  using  the  General 
Materials  (GEM)  diffractometer  at  the  ISIS  high-intensity 
pulsed  neutron  source  (Rutherford  Appleton  Laboratory, 
Didcot,  U.K.)[30]  and  yielded  the  neutron  total¬ 
scattering  structure  factors  (SSF),  S(Q): 


s(Q) 


1  da 
NdQ 


[2] 


Here  /  is  the  differential  neutron  cross-section  per 
unit  solid  angle  Q  for  the  3-component  alloy,  c,-  and  bt 
are,  respectively,  the  atomic  fraction  and  the  coherent 
bound  neutron  scattering  length  of  element  i.  Each  S(Q) 
was  Fourier  transformed  to  give  a  real  space  neutron 
total  RDF,  G(r),  defined  as[31,32] 


G{r)  =  Y  ciCjbibj[gij{r)  -  1] 
ij=  1 

oo 

=  2^-J  QS(Q)L(Q)MrQ)dQ  [3] 
0 


The  Lorch  modification  function/331  L(Q),  with  a 
maximum  momentum  transfer  Qmax  of  25  A  1  was  used 


III.  RESULTS  AND  DISCUSSION 

A.  Structure  Factors  and  Radial  Distribution  Functions 

The  experimental  SSFs,  S(Q),  for  the  six  amorphous 
alloys  are  given  in  Figure  1.  The  S(Q)  curves  show  a 
pronounced  peak  in  the  Q  range  from  2  to  3  A  1  and  several 
small  diffuse  peaks  at  higher  Q  values.  The  intensities  of 
these  peaks  decrease  quickly  with  Q  and  no  oscillations  are 
apparent  beyond  -15  A  /  As  the  Cu  and  Zn  concentra¬ 
tions  increase,  the  first  diffraction  peak  broadens  (Figure  1 ). 
The  presence  of  a  pre-peak  at  Q  «  0.8  to  1.5  A  1  may 
indicate  the  presence  of  MRO  in  these  alloys/34,351  as  well  as 
the  fluctuation  in  the  atomic  scattering  cross-sections:  the 
solutes,  Cu,  Zn  and  Mg,  have  larger  neutron  scattering 
cross-sections  than  the  solvent,  Ca. 

The  experimental  RDFs,  G(r),  are  given  in  Figure  2  as 
solid  lines.  Systematic  changes  of  the  shape  of  the  first 
RDF  peak,  which  is  located  at  r  «  2.1  to  4.5  A,  with  the 
alloy  composition  are  observed  (Figure  2).  For  example, 
the  first  RDF  peak  in  Ca60Mg15Cu25  has  a  maximum  at 
r  =  3.05  A  and  two  shoulders,  one  at  r  =  2.58  A  and 
another  at  r  =  3.70  A.  As  the  amount  of  Ca  decreases 
and  the  amounts  of  Cu  and  Mg  increase,  the  intensity  of 
the  first  shoulder  increases  and  the  intensity  of  the 
second  shoulder  decreases.  On  the  other  hand,  an 
increase  in  the  concentration  of  Mg  at  the  expense  of 
Cu  results  in  the  intensity  of  the  second  shoulder  to 
exceed  the  intensity  of  the  former  maximum  (compare 
RDFs  of  Ca60Mgi5Cu25  and  Ca60Mg25Cu15  in 
Figure  2).  Replacing  Cu  with  Zn  noticeably  decreases 
the  intensity  of  the  first  shoulder.  Because  the  first  RDF 
peak  is  a  sum  of  the  concentration-  and  scattering 
length-weighted  peaks  of  the  six  relevant  PRDFs 
(Eq.  [3]),  such  a  strong  concentration  dependence  of 
the  shape  of  the  first  RDF  peak  is  apparently  due  to 
different  partial  coordination  numbers  and  pair  bond 
distances  in  the  studied  alloys.  The  very  different 
neutron  scattering  length  of  Cu  ( bcu  —  7.718  fm)  and 
Zn  ( bZn  =  5.680  fm)  need  also  to  be  taken  into  account 
when  comparing  Ca-Mg-Cu  and  Ca-Mg-Zn. 
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Fig.  1  Total  neutron  scattering  factors  of  several  Ca  Mg  Zn  and 
Ca  Mg  Cu  amorphous  alloys. 

Figure  3  shows  six  QMD-simulated  PRDFs  for  the 
CawlMg15Zn25,  and  Ca6oMgi5Cu25  amorphous  alloys. 
The  PRDFs  are  rather  smooth  since  they  account  for 
atom  vibrations.  The  first  peak  of  each  of  the  PRDFs 
corresponds  to  the  atom  pair  interactions  within  the  first 
coordination  shell  and  it  provides  information  on  the 
nearest  neighbor  atom  pair  bond  distance  distributions 
and  partial  coordination  numbers  around  the  alloying 
elements.  Comparison  of  the  PRDFs  of  the 
Ca60Mg15Zn25,  and  Ca60Mg15Cu25  amorphous  alloys 
shows  that  the  intensities  of  the  Ca-Ca,  Ca-Mg  and 
Ca-Zn  peaks  in  the  CaeoMgisZ^s  alloy  are  similar  to 
the  intensities  of  the  respective  Ca-Ca,  Ca-Mg  and 
Ca-Cu  peaks  in  the  Ca«)Mgi5CU2s  alloy.  At  the  same 
time  the  intensities  of  the  Zn-Zn  and  Mg-Zn  peaks  are 
smaller  and  the  intensity  of  the  Mg-Mg  peak  is  larger  in 
Ca6oMgi5Zn25  than  the  respective  intensities  of  the 
Cu-Cu,  Mg-Cu  and  Mg-Mg  peaks  in  Ca6oMgi5Cu25. 
This  interesting  observation  may  indicate  that  the 
Ca-Ca,  Ca-Mg  and  Ca-Zn(Cu)  interactions  are  similar 
in  these  alloys,  while  the  tendencies  for  the  Cu-Cu  and 
Mg-Cu  bonding  are  stronger  than  for  the  Zn-Zn  and 
Mg-Zn  bonding.  Mg  atoms  seem  to  have  higher 
tendency  to  cluster  in  the  Ca60Mg15Zn25  alloy  than  in 
the  Ca6oMgi5Cu25  alloy. 

Figure  4  illustrates  the  composition  dependence  of 
two  PRDFs,  gcu-CuM  and  gca-Mg(r)  for  the  four  Ca-Mg- 
Cu  amorphous  alloys.  The  effect  of  the  composition  on 
the  position  of  the  first  PRDF  peak  is  weak  relative  to 
the  effect  on  the  positions  of  the  higher  order  peaks. 
This  indicates,  to  the  first  order,  that  bond  distances 
between  respective  atom  pairs  in  the  first  coordination 
shell  are  weakly  affected  by  the  composition  in  the 
Ca-Mg-Cu  system.  The  noticeable  composition  depen¬ 
dence  in  the  positions  of  the  second  and  higher  order 
PRDF  peaks  should  then  indicate  that  the  pair  bond 
angular  distributions  are  different  in  these  alloys.1'31 
A  stronger  composition  dependence  of  the  intensity  of 


Fig.  2  Experimental  (solid  lines)  and  QMD  simulated  (dashed 
lines)  total  radial  distribution  functions,  G(r),  of  Ca  Mg  Zn  and 
Ca  Mg  Cu  amorphous  alloys. 

the  first  peak  for  the  Cu-Cu  PRDF  than  for  the  Ca-Mg 
PRDF  may  indicate  that  the  Cu-Cu  chemical  order  is 
stronger  than  the  Ca-Mg  chemical  order. 

Using  the  QMD-simulated  PRDFs,  gj(r),  known 
concentrations,  c„  and  known  neutron  scattering 
lengths,  bj,  of  the  alloying  elements,1361  the  simulated 
total  RDFs  were  obtained  for  every  alloy  by  utilizing 
Eq.  [3]  and  compared  with  the  respective  experimental 
neutron  RDFs  (Figure  2).  A  very  good  agreement  is 
seen.  Since  the  QMD  simulations  were  conducted 
completely  independently  from  the  neutron  diffraction 
experiments,  these  excellent  concords  indicate  that  the 
simulated  PRDFs  are  statistically  consistent  with  the 
amorphous  structures  of  the  studied  alloys  and  they  can 
be  used  to  obtain  statistical  information  on  the  effect  of 
the  alloy  composition  on  the  SRO  in  these  alloy  systems. 

B.  Atom  Pair  Bond  Distances 

Using  the  QMD-simulated  PRDFs,  the  most  frequent 
(mode,  r¥ode),  weighted  average  (mean,  /-¥ean)  and  cut¬ 
off  (maximum,  r 0“  off)  distances  between  atom  pairs  in 
the  first  coordination  shell  were  calculated  for  the 
studied  alloys  and  their  values  are  given  in  Tables  II, 
III,  and  IV,  respectively.  The  weighted  average  dis¬ 
tances,  /-¥eilI\  between  the  i  and  j  atoms  were  calculated 
using  the  following  equation12  '1: 

jCut  off  -Cut  off 

u  r# 

rj}lean  =  J  r<jgij(r)dr/  J  gtj{r)dr  [5] 
o  o 

The  /-Mode  values  (Table  II)  are  always  smaller  than 
the  respective  rMcan  values  (Table  III).  This  indicates 
that  the  first-shell  gi,{r)  peaks  are  asymmetric  (non- 
Gaussian).  The  difference  is  larger  for  the  Zn-  than  for 
the  Cu-containing  pairs,  which  may  indicate  a  stronger 
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(a) 

Fig.  3  QMD  simulated  partial  radial  distribution  functions,  gg(r),  for 


(a) 

Fig.  4  Composition  dependence  of  two  PRDFs:  (a)  go,  o, (r)  and  ( b ) 

interaction  of  Cu,  than  Zn,  with  other  elements. 
The  rM<Kic  and  rMean  values  can  be  compared  with 
atomic  separations  in  crystalline  metals  and  alloys. 
Crystalline  metallic  separations,  rM,  and  covalent  dis¬ 
tances,  rc,  were  estimated  as  a  sum  of  metallic  and 
covalent  radii  of  respective  elements  (Table  V).  The 
metallic  radii  are  half  the  nearest  distance  between 
atoms  in  crystal  lattices  of  the  respective  pure  metals.[37] 
The  covalent  atomic  separations  were  obtained  from 
Reference  38,  based  on  the  analysis  of  the  nearest 
distances  in  crystalline  intermetallic  compounds  re¬ 
ported  in  the  Cambridge  Structural  Database.  From 
T able  V,  the  rMcan  and  /‘Mode  values  for  Ca-Cu,  Ca-Zn, 
and  Mg-Cu  atom  pairs  are  smaller  than  rM.  However, 
for  all  other  atomic  pairs  /"Mean  values  are  larger  and 
'"Mode  values  are  smaller  than  r\ j.  At  the  same  time, 
'"Mode  f°r  the  atom  pairs  containing  Cu  are  even  smaller 
than  the  respective  covalent  bond  distances  in  crystalline 
compounds.  The  Mg-Cu  (2.73  A),  Ca-Cu  (3.03  A)  and 
Ca-Zn  (3.14  A)  mode  pair  distances  are  much  shorter 
than  those  calculated  from  the  Ca-Ca,  Mg-Mg,  Cu-Cu 
and  Zn-Zn  rModc  values  (i.e.,  2.82,  3.13,  and  3.20  A, 
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(b) 

Ca60Mg15Zn25  and  (b)  Ca6oMg|5Cu25  amorphous  alloys. 


(b) 

Mg (r)  for  Ca  Mg  Cu  amorphous  alloys. 

respectively,  Table  V).  Shortening  of  Mg-Cu  and  Y-Cu 
bonds  has  also  been  observed  in  amorphous 
MgfioCusoYm-1391  Such  bond  shortening  often  indicates 
strong  chemical  interactions  between  these  atom  pairs. 
No  bond  shortening  is  observed  between  Ca-Mg  and 
Mg-Zn  pairs. 

The  observation  that  some  nearest-neighbor  bond 
lengths  are  shorter  than  those  in  the  competing  long- 
range  ordered  crystals  may  help  explain  why  these 
Ca-Mg-Cu  and  Ca-Mg-Zn  alloys  are  easy  glass  formers. 
The  equilibrium  crystalline  phases  that  compete  with  the 
amorphous  phase  appear  to  minimize  total  energy  by 
sacrificing  optimal  SRO  to  achieve  beneficial  long-range 
order.  On  the  other  hand,  the  absence  of  long-range 
atomic  order  allows  metallic  glasses  to  reduce  the  energy 
difference  between  the  metastable  glass  and  equilibrium 
crystal(s)  by  optimizing  short-range  atomic  interactions 
and  arranging  the  atoms  in  efficiently  packed  clus¬ 
ters.140  43)  The  largest  contribution  to  condensed  phase 
stability  is  expected  to  come  from  nearest-neighbor 
interactions  and  volume  minimization,  so  that  the  opti¬ 
mized  short-range  interactions,  as  well  as  the  presence  of 
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MRO,  in  metallic  glasses  can  produce  structures  that 
have  only  a  small  energetic  disadvantage  relative  to  the 
equilibrium  crystalline  structure. [22  2T|  Kinetic  con¬ 
straints  from  quenching  and  packing  frustration  restrict 
the  long  range  atomic  redistribution  needed  to  achieve 
long-range  order  and  to  further  minimize  the  total 


Table  II.  The  Most  Frequent  (Mode)  Distances  (in  Angstroms) 
Between  Pair  Atoms  in  Ca-Mg-TM  (TM  =  Cu  or  Zn) 
Amorphous  Alloys,  in  Accord  to  QMD  Simulations 


Alloy 

Ca 

Ca 

Ca 

Mg 

Ca 

TM 

Mg 

Mg 

Mg 

TM 

TM 

TM 

CaSoMgi5Zn25 

3.81 

3.44 

3.14 

3.16 

2.83 

2.67 

Ca60Mg25Ziii5 

3.82 

3.48 

3.14 

3.18 

2.94 

2.57 

Ca60Mgi5Cu25 

3.78 

3.43 

3 

3.14 

2.72 

2.54 

CaS0Mg25Cu15 

3.78 

3.46 

3.03 

3.08 

2.72 

2.51 

Ca50Mg25Cu25 

3.75 

3.41 

3.07 

3.13 

2.75 

2.48 

Ca40Mg25CU35 

3.75 

3.42 

3.02 

3.14 

2.71 

2.47 

Table  III.  Weighted  Average  (Mean)  Distances 
(in  Angstroms)  Between  Pair  Atoms  in  the  First  Shell 
in  Ca-Mg-TM  (TM  =  Cu  or  Zn)  Amorphous  Alloys, 
in  Accord  to  QMD  Simulations 


Alloy 

Ca 

Ca 

Ca 

Mg 

Ca 

TM 

Mg 

Mg 

Mg 

TM 

TM 

TM 

Ca60Mgi5Zn25 

4.00 

3.67 

3.30 

3.27 

3.06 

2.83 

Ca60Mg25Ziii5 

4.01 

3.65 

3.31 

3.37 

3.05 

2.76 

CaSoMgi5Cu25 

4.00 

3.67 

3.18 

3.25 

2.88 

2.60 

CaS0Mg25Cu15 

4.00 

3.64 

3.19 

3.29 

2.87 

2.60 

Ca5oMg25Cu25 

4.01 

3.62 

3.25 

3.24 

2.89 

2.58 

Ca40Mg25CU35 

3.95 

3.64 

3.21 

3.31 

2.83 

2.59 

Table  IV.  The  First  Shell  Cut-Off  Distances  Between  Pair 
Atoms  (in  Angstroms)  in  Ca-Mg-TM  Amorphous  Alloys 
(TM  =  Cu  or  Zn),  in  Accord  to  QMD  Simulations 


Alloy 

Ca 

Ca 

Ca 

Mg 

Ca 

TM 

Mg 

Mg 

Mg 

TM 

TM 

TM 

CaSoMgi5Zn25 

4.91 

4.55 

4.21 

4.33 

4.13 

3.46 

CaS0Mg25Zni5 

4.95 

4.55 

4.2 

4.33 

3.85 

3.46 

Ca60Mgi5Cu25 

4.93 

4.72 

4.19 

4.01 

3.84 

3.33 

Ca60Mg25Cui5 

4.94 

4.58 

4.13 

4.05 

3.89 

3.31 

Ca50Mg25Cu25 

4.99 

4.58 

4.4 

3.97 

3.84 

3.42 

Ca40Mg25Cu35 

4.88 

4.63 

4.26 

4.22 

3.58 

3.25 

system  energy,  thus  favoring  glass  formation.  The 
present  results  suggest  that  metallic  glasses  not  only 
have  SRO,  but  may  have  "better”  combination  of  short- 
range  topological  and  chemical  order  (in  terms  of 
optimal  bond  length  and/or  atomic  arrangement)  than 
the  competing  crystals. 


C.  Coordination  Numbers  and  Chemical  Short 
Range  Order 

The  average  partial  and  total  coordination  numbers 
around  Ca,  Mg  and  TM  atoms  (TM  is  Cu  or  Zn),  from  the 
QMD  simulations,  are  given  in  Table  VI.  There  is  almost 
no  effect  from  replacing  Cu  with  Zn  on  the  partial  and  total 
coordination  numbers  around  Ca  and  Mg  atoms  in  the 
Ca60Mg15TM25  and  Ca60Mg25TM15  amorphous  alloys 
(the  first  4  columns  in  Table  VI).  The  total  coordination 
number  around  Ca  and  Mg  are  CNCa  =  14.3  and 
CNMg  =  11.3  to  11.5  for  Ca60Mg25TM15  and  CNCa  = 

14.4  and  CNMg  =  11.9  to  12.0  for  Ca60Mg15TM25.  The 
number  of  TM  and  Ca  increase  and  Mg  decrease  in  the  first 
coordination  shell  of  Ca  and  Mg  clusters  with  an  increase 
in  the  concentration  of  the  TM  element  in  these  alloys.  At 
the  same  time,  Cu  has  a  slightly  different  neighbor 
environment  than  Zn.  In  particular,  the  amount  of  Ca  is 
slightly  lower  and  Mg  slightly  higher  around  Cu  than 
around  Zn.  This  is  in  agreement  with  the  earlier  discussed 
observation  that  the  Mg-Cu  interactions  are  stronger  than 
Mg-Zn  interactions. 

Increasing  the  Cu  concentration  from  15  to  35  at.  pet 
continuously  increases  the  total  coordination  numbers 
CNCa  from  14.3  to  16.2  and  CNMg  from  11.3  to  12.8, 
due  to  a  faster  increase  in  the  number  of  smaller  Cu 
atoms  (CNCa_Cu  and  CNMg_Cu  increase  from  1.6  to  4.9 
and  from  1.1  to  3.2,  respectively)  and  a  slower  decrease 
in  the  amount  of  the  larger  Ca  atoms  (CNCa_Ca 
decreases  from  9.2  to  7.0  and  CNMg_Ca  decreases  from 
8.5  to  6.9)  (Table  VI).  On  the  other  hand,  CNcu  almost 
does  not  depend  on  the  alloy  composition.  However, 
similar  to  the  environment  around  Ca  and  Mg,  CNCu_Cu 
increases  from  0.6  to  1.6  and  CNCu-Ca  decreases  from 

6.4  to  5.6  with  the  increase  in  Cu  concentration.  The 
number  of  Mg  atoms  around  Ca,  Mg  and  Cu  has  a 
tendency  to  weakly  increase  with  an  increase  in  the 
concentration  of  Cu  (Table  VI). 

To  determine  if  the  atoms  in  the  amorphous  structure 
are  distributed  randomly  or  they  have  a  tendency  to 
form  specific  clusters,  the  degree  of  chemical  short  range 


Table  V.  Metallic,  rm,|3?l  Covalent,  i'c,[3S|  mode,  rMode»  and  mean,  i'Mean,  Bond  Distances  Between  Ca,  Mg  and  Cu(Zn)  Atom  Pairs 


Ca  Ca 

Ca  Mg 

Ca  Cu 

Mg  Mg 

Mg  Cu 

Cu  Cu 

Ca  Zn 

Mg  Zn 

Zn  Zn 

i'm  (A) 

3.94 

3.57 

3.25 

3.20 

2.88 

2.56 

3.31 

2.94 

2.68 

rc  (A) 

3.52 

3.17 

3.08 

2.82 

2.73 

2.64 

2.98 

2.63 

2.44 

^  Mode  (A) 

3.78 

3.44 

3.03 

3.14 

2.73 

2.50 

3.14 

2.89 

2.62 

^  Mean  (A) 

4.00 

3.65 

3.21 

3.29 

2.87 

2.59 

3.31 

3.06 

2.80 

100  pet  (/*Mean/^Mode  1) 

5.82 

6.10 

5.94 

4.78 

5.13 

3.60 

5.41 

5.88 

6.87 

100  pet  (rM/l*Mode  1) 

4.23 

3.78 

7.26 

1.91 

5.49 

2.40 

5.41 

1.73 

2.29 

100  pet  (rc/rM ode  1) 

6.88 

7.85 

1.65 

10.19 

0.00 

5.60 

5.10 

9.00 

6.87 

100  pet  (/'M/rMean  1) 

1.50 

2.19 

1.25 

2.74 

0.35 

1.16 

0.00 

3.92 

4.29 

rm  and  rc  are  compared  to 

'Mode  and  '"Mean  from  the  QMD  PRDFs. 

■  ''Mode  and  rMean  are  the  average  values  for  the  studied  alloys. 
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Table  VI.  Partial  and  Total  Coordination  Numbers  Around  Ca,  Mg  and  TM  (TM  =  Cu  or  Zn)  Atoms  in  the  Ca-Mg-Zn 
and  Ca-Mg-Cu  Amorphous  Alloys,  in  Accord  to  QMD  Simulations 


CN 

Ca6oMgi5Zn25 

Ca60Mg25Zni5 

Ca60Mg15Cu25 

Ca6oMg25Cu15 

Ca50Mg25CU25 

Ca40Mg25Cu35 

Ca  Ca 

9.3 

9.1 

9.3 

9.2 

8.6 

7.0 

Ca  Mg 

2.1 

3.4 

2.3 

3.5 

3.8 

4.3 

Ca  TM 

3.0 

1.8 

2.8 

1.6 

3.2 

4.9 

Total 

14.4 

14.3 

14.4 

14.3 

15.6 

16.2 

Mg  Ca 

8.5 

8.3 

9.0 

8.5 

7.7 

6.9 

Mg  Mg 

1.5 

2.4 

1.0 

1.8 

2.1 

2.6 

Mg  TM 

1.9 

0.8 

2.0 

1.1 

2.1 

3.2 

Total 

11.9 

11.5 

12.0 

11.3 

11.8 

12.8 

TM  Ca 

7.3 

7.1 

6.7 

6.4 

6.4 

5.6 

TM  Mg 

1.1 

1.4 

1.2 

1.8 

2.1 

2.3 

TM  TM 

0.8 

0.8 

1.0 

0.6 

1.0 

1.6 

Total 

9.2 

9.3 

8.9 

8.8 

9.5 

9.4 

Table  VII.  The  Short-Range  Order  Parameters,  a,w, 
for  the  First  Coordination  Shell  of  Ca,  Mg  and  TM  Atoms 
in  the  Ca-Mg-TM  Amorphous  Alloys  (TM  is  Cu  or  Zn) 


ai(Jk)'  i 

Ca 

Mg 

TM 

Ca6oMgi5Zn25 

0.11 

0.03 

0.22 

Ca6oMg25Zn15 

0.09 

0.06 

0.08 

Ca6oMgi5Cu25 

0.11 

0.08 

0.18 

Ca6oMg25Cu15 

0.11 

0.13 

0.10 

Ca5oMg25Cu25 

0.10 

0.11 

0.19 

Ca40Mg25CU35 

0.05 

0.05 

0.29 

order  (CSRO)  in  the  first  coordination  shell  of  Ca,  Mg 
and  TM  atoms  was  estimated  using  a  CSRO  param¬ 
eter,  oiiQk),  which  was  defined  for  ternary  alloys  in 
Reference1441: 

a,-(/fc)  =  1  -  (CNy  +  CNjfc) / ( (cj  +  ck) CNj )  for  ifjfk 

[6] 

Here  Cj  is  the  atomic  fraction  of  the  element  j,  C  N  ,y  is 
the  partial  coordination  number  of  element  j  in  the  first 
coordination  shell  of  the  element  /,  and  CN,-  is  the  total 
coordination  number  of  element  /.  Negative  values 
indicate  the  presence  of  CSRO,  i.e.,  the  total  concen¬ 
tration  of  j  and  k  atoms  in  the  first  coordination  shell 
exceed  their  total  average  concentration  in  the  alloy. 
The  positive  values  indicate  the  excess  of  /-elements  in 
the  first  coordination  shell  of  /  atom,  i.e.,  chemical  short 
range  clustering  (CSRC)  of  atoms  /. 

The  y.Kjk)  values  for  the  studied  Ca-Mg-TM  glasses  are 
given  in  Table  VII.  The  parameter  a,y/c)  is  positive  for 
Ca  atoms  and  is  negative  for  Mg  and  TM  atoms.  This 
result  indicates  that  Ca  atoms  tend  to  cluster  with  each 
other,  while  Mg  and  TM  atoms  have  a  tendency  to 
attract  unlike  atoms  in  the  first  coordination  shell  and 
form  CSRO.  This  tendency  for  CSRO  is  stronger  for  the 
TM  elements  than  for  Mg  in  all  studied  alloys,  except 
Ca6oMg25Cu15.  The  CSRO  of  the  TM  element  with  Ca 
and  Mg  increases  with  an  increase  in  the  concentration 
of  the  TM  element  in  the  alloy.  The  presence  of  CSRO 
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near  TM  and  Mg  and  CSRC  near  Ca  atoms  indicates 
that  the  amorphous  structure  of  the  Ca-Mg-TM  alloys  is 
not  a  random  solution  of  the  alloying  elements,  but  it 
consists  of  TM-  and  Mg-centered  atomic  clusters,  with 
specific  CSROs.  Stronger  CSRO  around  the  TM  ele¬ 
ments  may  indicate  that  the  TM-centered  clusters 
dominate  over  the  Mg-centered  clusters. 

D.  Voronoi  Analysis  of  the  Atomic  Clusters 
in  Ca-Mg-TM  Amorphous  Alloys 

Voronoi  tessellation  and  the  resulting  Voronoi  cells 
allow  identification  of  characteristic  coordination  poly- 
hedra  present  in  an  amorphous  structure.  23,421  A  coor¬ 
dination  polyhedron  is  defined  as  an  /-centered  cluster 
with  vertices  at  the  first-shell  atom  positions  and  edges 
coinciding  with  the  interatomic  bonds  in  the  first 
shell.1451  Each  coordination  polyhedron  can  be  assigned 
a  Voronoi  signature  (n3,  «4,  n5,  n6),  where  nm  is  the 
number  of  vertices  common  to  m  polyhedron  faces  (or 
edges).1461  m  is  also  called  the  vertex  coordination.  Those 
/-centered  clusters  with  the  same  Voronoi  signature  are 
considered  to  be  topologically  equivalent  (even  though 
they  may  not  be  identical). 

Table  VIII  lists  Ca,  Mg  and  TM-centered  clusters 
identified  in  the  amorphous  structures  of  the  studied 
Ca-Mg-TM  alloys.  Only  clusters,  which  are  present  at  a 
fraction,  /P,  of  0. 1  or  higher  in  at  least  one  of  the  alloys  are 
tabulated.  Just  three  types  of  the  Ca-centered  polyhedra 
have  the  fractions /P  >  0.1.  These  are  (0,2, 8, 4),  (0,1,10,4) 
and  (0,2, 8, 5).  The  total  fraction  of  these  clusters  varies 
from  0.36  in  Ca60Mg25Zn15  to  0.075  in  Ca40Mg25Cu35  and 
many  other  different  types  of  Ca-centered  polyhedra  are 
present  at  /P  -  0.05  or  less.  Many  different  types  of 
clusters,  with  low  fraction  of  each  of  the  type,  may 
indicate  that  the  Ca-centered  clusters  are  not  likely 
building  blocks  for  the  amorphous  structures. 

Six  characteristic  Mg-centered  and  five  TM-centered 
clusters  with  /P  >  0.1  are  also  shown  in  Table  VIII.  All 
of  them  except  (1,3, 3, 3)  are  Kasper-type  polyhedra,1451 
which  account  for  36  to  80  pet  among  the  Mg-centered 
clusters  and  44  to  82  pet  among  the  Cu-centered  clusters 
(these  clusters  are  the  least  populated  in  Ca40Mg25Cu35 
and  most  populated  in  Ca60Mg15Zn25).  The  high  fraction 

VOLUME  44A,  MAY  2013  1985 


6 

Distribution  Statement  A.  Approved  for  public  release;  distribution  unlimited. 


Table  VIII.  Fraction  of  Clusters  Identified  in  the  Amorphous  Structures  of  Ca-Mg-TM  Alloys 


Cluster  type 

Ca6oMg15Zn25 

Ca60Mg25Zn15 

Ca60Mg15Cu25 

Ca60Mg25Cu15 

Ca50Mg25Cu25 

Ca40Mg25Cu35 

Ca  centered  clusters 

(0,2, 8, 4) 

0.12 

0.15 

0.08 

0.10 

0.05 

0.01 

(0,1,10,4) 

0.04 

0.11 

0.18 

0.11 

0.07 

0.04 

(0,2,8, 5) 

0.04 

0.10 

0.06 

0.04 

0.10 

0.03 

Mg  centered  clusters 

(0,2,8, 1) 

0.23 

0.10 

0.13 

0.16 

0.06 

0.00 

(0, 3,6,3) 

0.10 

0.04 

0.13 

0.10 

0.02 

0.00 

(0,0,12,0) 

0.20 

0.14 

0.07 

0.04 

0.06 

0.08 

(0,2, 8,2) 

0.10 

0.18 

0.03 

0.10 

0.10 

0.02 

(0, 3,6,4) 

0.07 

0.16 

0.17 

0.00 

0.16 

0.06 

(0,1,10,2) 

0.10 

0.06 

0.03 

0.08 

0.04 

0.20 

TM  centered  clusters 

(0, 3,6,0) 

0.30 

0.13 

0.24 

0.27 

0.16 

0.11 

(1, 3,3,3) 

0.08 

0.00 

0.10 

0.02 

0.06 

0.07 

(0, 3,6,1) 

0.20 

0.10 

0.14 

0.10 

0.10 

0.04 

(0, 2,8,0) 

0.16 

0.27 

0.12 

0.03 

0.12 

0.11 

(0,2,8, 1) 

0.08 

0.20 

0.06 

0.03 

0.04 

0.10 

Only  clusters  populated  with  the  fractions  of  0.1  or  higher  in  at  least  one  of  the  alloys  are  tabulated. 

Table  IX.  Global  and  Atomic  Packing  Fractions  in  Ca-Mg-TM  Amorphous  Alloys 

Alloy 

Packing  Fraction 

Average  Voronoi  Volume  (A3) 

Atomic  Packing  Fraction 

Ca 

Mg 

TM 

Ca 

Mg 

TM 

Ca60Mgi5Zn25 

0.76 

41.0 

24.2 

14.3 

0.78 

0.71 

0.71 

Ca60Mg25Ziii5 

0.76 

40.6 

25.3 

15.9 

0.79 

0.68 

0.64 

Ca60Mg15Cu25 

0.78 

39.9 

23.7 

12.4 

0.80 

0.72 

0.71 

CaSoMg25Cu15 

0.77 

40.5 

24.2 

12.8 

0.79 

0.71 

0.68 

Ca50Mg25Cu25 

0.79 

38.8 

24.1 

12.6 

0.83 

0.71 

0.70 

Ca40Mg25CU35 

0.80 

38.0 

22.6 

12.0 

0.84 

0.76 

0.73 

of  the  selected  TM-  and  Mg-centered  Kasper-type 
polyhedra  should  indicate  that  the  amorphous  structure 
of  the  Ca-Mg-TM  alloys  mainly  consists  of  these 
clusters,  which  is  in  agreement  with  the  recently  pro¬ 
posed  ECP  model  of  the  amorphous  structure.  ,47] 
Since  Kasper  polyhedra  are  polytetrahedral,  and  many 
other  polyhedra  also  contain  tetrahedra,  we  conclude 
that  the  topological  SRO  of  Ca-Mg-TM  is  polytetrahe¬ 
dral  in  nature,  which  has  been  demonstrated  in  many 
other  metallic  glasses.1221  The  current  study  therefore 
further  supports  the  idea  that  the  SRO  of  MGs  is 
characterized  by  polytetrahedral  packing  via  Kasper 
polyhedra  and  their  distorted  variants. 

The  analysis  of  the  most  populated  clusters  given  in 
Table  VIII  indicates  that  five-coordinated  vertices 
dominate  in  all  the  clusters;  however,  the  amorphous 
structures  of  the  Ca-Mg-TM  glasses  contain  a  low 
fraction  of  icosahedral  SRO  (Kasper  polyhedra  with 
(0,0,12,0)  Voronoi  index).  This  is  presumably  due  to 
the  large  mismatch  of  atomic  sizes  in  the  nearest 
neighbor  shell  (e.g.,  rCa  =  1.97  A  and  rCu  =  1.28  A), 
which  renders  the  distorted  Kasper  polyhedra  more 
topologically  favored.  A  large  fraction  of  five-fold 
bonds  and  the  lack  of  icosahedral  SRO  have  also  been 
found  in  the  MD-simulated  amorphous  structure  of 
Mg-Cu  alloys. [24]  The  dominance  of  five-fold  vertices 
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suggests  that  tetrahedra  prefer  to  cluster  into  pentag¬ 
onal  bi-pyramids. 

Table  VIII  also  demonstrates  the  different  local 
topology  in  Ca6oMg15Zn25  and  Ca60Mg15Cu25  (columns 
1  and  3),  which  differs  by  only  the  type  of  the  TM 
element.  Specifically,  because  Cu  is  smaller  in  size,  and 
has  obvious  bond  shortening  with  both  Ca  and  Mg,  it  is 
expected  to  induce  higher  CN  around  Ca/Mg,  as  well  as 
more  distortion  in  the  Kasper  polyhedra.  Indeed,  the 
most  popular  cluster  type  around  Ca  is  (0,1,10,4)  with 
CN  =  15  in  Ca60Mg15Cu25,  and  (0,2, 8,4)  with  CN  =  14 
in  Ca60Mgi5Zn25.  For  Mg  centered  clusters,  it  is  (0,3, 6, 4) 
with  CN  =  13  in  Ca60Mg15Cu25,  and  (0,2,8, 1)  with 
CN  =  11  and  (0,0,12,0)  with  CN  =  12  in  Ca60Mg15Zn25. 

E.  Local  Volume  and  Packing  Fraction 

Assuming  atoms  are  hard  spheres,  the  global  packing 
fraction  can  be  easily  calculated  using  the  atomic  radii 
and  the  density  of  the  alloy.  This  approach  yields  a 
packing  fraction  of  0.76  to  0.80  for  the  alloys  in  this 
study  (Table  IX).  However,  this  treatment  is  very  crude, 
and  it  does  not  provide  information  on  local  packing 
fraction  specific  to  an  atom,  an  element,  or  a  cluster.  In 
the  following  we  propose  a  method  that  can  be  used  to 
describe  the  local  volume  and  packing  fraction. 
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Table  X.  Cluster  Packing  Fractions  in  Ca-Mg-TM  Amorphous  Alloys 

Average  Vertex  Coordination  Average  Cluster  Volume  (A3)  Cluster  Packing  Fraction 


Alloy 

Ca 

Mg 

TM 

Ca 

Mg 

TM 

Ca 

Mg 

TM 

Ca60Mg15Zn25 

5.13 

4.97 

4.78 

107.6 

83.3 

64.0 

0.77 

0.75 

0.76 

Ca6oMg25Zn15 

5.14 

5.00 

4.84 

108.5 

84.0 

65.1 

0.77 

0.74 

0.74 

Ca60Mg15Cu25 

5.14 

5.00 

4.77 

104.2 

81.9 

57.9 

0.79 

0.77 

0.77 

Ca6oMg25Cu15 

5.12 

4.96 

4.77 

107.7 

82.0 

58.7 

0.78 

0.76 

0.76 

Ca5oMg25Cu25 

5.17 

5.02 

4.80 

103.7 

77.6 

58.4 

0.80 

0.77 

0.78 

Ca4oMg25Cu35 

5.16 

5.04 

4.86 

96.9 

73.6 

53.4 

0.82 

0.80 

0.80 

The  weighted  Voronoi  tessellation  assigns  each  atom 
with  a  Voronoi  cell,  and  its  volume  can  be  considered  as 
the  atomic  volume  (i.e.,  the  volume  occupied  by  the 
center  atom  at  its  position).  The  summation  of  the 
atomic  volume  defined  as  such  equals  to  the  total 
volume  of  the  alloy.  Therefore,  the  ratio  between  the 
hard  sphere  volume  (calculated  from  the  metallic  radius) 
and  the  atomic  volume  (Voronoi  cell  volume)  is  on 
average  smaller  than  1,  and  reflects  the  packing  density 
around  the  specific  atom  (i.e.,  atomic  packing  fraction  in 
Table  IX).  The  challenge  now  is  how  to  define  the 
packing  fraction  for  a  short-range  cluster,  and  the  key  is 
to  find  out  a  way  to  quantify  the  actual  volume  occupied 
by  a  cluster  in  the  amorphous  matrix.  An  /-centered 
cluster  consists  of  the  atom  /  in  the  center  and 
CN,-  =  Ylj=i  CN y  atoms  in  the  first  shell.  The  atoms  in 
the  first  shell  are  shared  with  other  neighbor  clusters. 
The  average  volume  of  the  cluster,  Vh  is  then  defined  as: 


CIV- 

vt  =  v/  +  2^ 


m1  +  k» 


[V] 


Here  v*  and  Vj  are  the  average  Voronoi  volumes  of 
atoms  i  and  j,  and  Ky  is  the  average  number  of  clusters, 
which  are  the  neighbors  to  the  /-centered  cluster  and 
which  share  the  atom  j  in  their  first  shell.  The  averaging 
is  conducted  over  the  whole  number  of  clusters  of  type  / 
presented  in  the  amorphous  structure  and  vertices  of 
type  j  in  the  cluster  /.  It  is  expected  that  K y  =  1  +  fiy, 
where  jiy  is  the  average  coordination  of  vertex  j  in  the 
/-centered  cluster: 


and  thus 


%  =  H  mufmi 

m=  3 


v>  ~v‘'  +  ^2 


i=  i 


CNyVj 
2  +  % 


[8] 

[9] 


Here  fmil  is  the  fraction  of  j  vertices  with  the  vertex 
coordination  my  in  all  /-centered  clusters.  Indeed,  in 
general  two  neighbor  clusters  can  share  a  face,  an  edge, 
or  a  vertex.  The  average  number  of  faces  or  edges,  which 
have  a  common  vertex  /  is  equal,  by  definition,  to  fiy.  A 
geometric  analysis  shows  that  the  neighbor  clusters 


having  a  common  vertex  j,  can  fill  space  only  if  they 
share  /-vertex-containing  faces  with  the  /-centered  clus¬ 
ter  and  thus  the  number  of  these  neighbors  should  be 
equal  to  pt,y.  One  additional  neighbor  cluster  will  share 
only  the  vertex  j  with  the  /-centered  cluster.  The  /re¬ 
values  can  be  directly  obtained  from  the  results  of  the 
Voronoi  tessellation  of  the  amorphous  structure. 

Equations  [8]  and  [9]  were  used  to  calculate  the  average 
vertex  coordinations  and  average  volumes  of  the  Ca-, 
Mg-  and  TM-centered  clusters  for  all  the  studied  alloys. 
The  results  are  given  in  Table  X.  The  average  coordina¬ 
tion  of  the  Ca-,  Mg-,  and  TM-vertices  is -5.1,  5.0,  and  4.8, 
respectively,  and  it  weakly  depends  on  the  alloy  compo¬ 
sition.  At  the  same  time,  the  average  cluster  volumes  are 
noticeably  composition-dependent.  As  it  is  expected  from 
the  numbers  and  radii  of  the  atoms  in  the  clusters,  the  Ca- 
centered  clusters  have  the  largest  volume  (-97  to  109  A3) 
and  the  TM-centered  clusters  have  the  smallest  volume 
(-53  to  64  A3).  The  Mg-centered  clusters  have  an 
intermediate  volume  of -74  to  84  A3.  The  cluster  volumes 
in  the  Zn-containing  alloys  are  larger  than  in  the 
respective  Cu-containing  alloys,  probably  because  the 
atomic  radius  of  Zn  is  larger  than  that  of  Cu,  and  Cu 
incurs  more  significant  bond  shortening.  In  the  Ca-Mg- 
Cu  alloys,  the  volume  of  the  clusters  decreases  with  an 
increase  in  the  concentration  of  Cu  and/or  decrease  in  the 
concentration  of  Ca/Mg  due  to  an  increasing  number  of 
Cu  and  decreasing  number  of  Ca/Mg  in  the  first  shell  of 
all  three  type  clusters. 

The  average  packing  fraction  of  the  /-type  clusters  was 
estimated  by  calculating  the  total  volume  of  the  atoms 
belonging  to  the  averaged  /-cluster  (using  the  equation 
similar  to  Eq.  [9],  where  the  Voronoi  volumes  are  replaced 
with  respective  atom  volumes)  and  dividing  it  by  Vfvet. 
Atomic  radii  reported  in  Table  V  were  used  to  calculate 
the  atom  volumes.  The  results  are  given  in  Table  X.  The 
Ca-centered  clusters  are  packed  slightly  more  densely 
than  the  Mg-  and  Cu-centered  clusters.  In  Ca-Mg-Zn 
alloys,  the  average  packing  fraction  of  Ca-centered 
clusters  is  -0.77  and  it  decreases  to  -0.74  to  0.75  for 
Mg-centered  and  to  -0.74  to  0.76  for  Zn-centered 
clusters.  The  packing  fraction  of  the  clusters  is  slightly 
smaller  in  the  Zn-containing  alloys  than  in  the  respective 
Cu-containing  alloys,  probably  due  to  stronger  interac¬ 
tions  of  Cu  with  other  elements  leading  to  more  intense 
bond  shortening.  An  increase  in  the  concentration  of  Cu 
in  the  Ca-Mg-Cu  alloys  increases  the  clusters’  packing 
fraction,  thus  supporting  the  presence  of  strong  Cu-Cu, 
Cu-Mg  and  Cu-Ca  bonding  in  these  amorphous  alloys. 
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Fig.  5  Site  and  orbital  projected  electronic  DOS  for  (a)  Ca60Mg15 
Zn25  and  ( b )  Ca60Mgl5Cu25  amorphous  alloys. 

The  cluster  volume  and  packing  fraction  data  do  not 
appear  to  be  systematically  correlated  with  the  GFA  of 
these  alloys.  This  is  not  surprising  because,  in  addition 
to  the  uncertainties  in  determining  the  volume  at  atomic 
level,  the  many  factors  that  can  influence  the  GFA, 
including  the  relative  stability  of  the  competing  crystal¬ 
line  phases  and  the  interfacial  energy,  often  make  it 
difficult  to  establish  a  simple  method  to  predict  the 
ability  of  an  alloy  to  form  glass. 

F.  Electronic  Structure  of  Ca-Mg-TM 
Amorphous  Alloys 

It  has  been  shown  in  the  above  discussion  that  the 
presence  of  Cu  or  Zn  in  the  Ca-Mg-TM  amorphous 
alloys  leads  to  different  bonding  characters  and  thus 
properties.  It  is  therefore  intriguing  to  examine  the 
electronic  structure  of  Ca-Mg-Cu  and  Ca-Mg-Zn,  espe¬ 
cially  the  differences  associated  with  the  TM  element.  To 
this  end,  site-  and  orbital-projected  electronic  density  of 
states  (DOS)  was  calculated  using  the  configurations 
obtained  in  QMD  simulation,  on  a  3x3x3 
Monkhorst-Pack  grid.  Only  valence  electrons  (2  for 
Ca,  2  for  Mg,  11  for  Cu,  and  12  for  Zn)  are  accounted 
for  as  the  low-lying  core  electrons  are  unlikely  to 
participate  in  the  bonding.  Wigner  radii  for  the  decom¬ 
position  are  1.947  A  for  Ca,  1.524  A  for  Mg,  1.270  A 
for  Zn,  and  1.312  A  for  Cu.  The  results  are  compared  in 
Figure  5. 

The  most  prominent  difference  between  Figure  5(a) 
(Ca6oMg15Zn25)  and  Figure  5(b)  (Ca60Mg15Cu25)  is  the 
very  different  position  of  the  narrow  and  intense 
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3d-band:  for  Zn  it  is  over  7  eV  below  the  Fermi  level, 
while  for  Cu  it  is  2  to  3  eV  below  the  Fermi  level.  As  a 
result,  the  3d-band  of  Zn  barely  overlaps  with  the  5-  and 
p-conduction/bonding  bands  in  Ca/Mg,  whereas  for  Cu 
the  overlap  is  significant.  Such  overlap  between  the 
3d-band  of  TM  and  the  sp-band  of  simple  metals  is 
known  to  result  in  strong  chemical  interactions  and 
orbital  hybridization,  as  previously  demonstrated  in 
Al-Cu,  where  the  Al-Cu  bonds  are  shortened  by  over 
6  pct.[20]  The  DOS  in  Figure  5  may  explain  why  Cu 
induces  more  obvious  bond  shortening  than  Zn  ( e.g ., 
Ca-Cu  bond  is  shortened  by  7.26  pet,  while  Ca-Zn  bond 
is  shortened  by  only  5.41  pet). 


IV.  CONCLUSIONS 

1.  The  atomic  structures  of  six  Ca-Mg-TM  ternary 
metallic  glasses  (TM  is  Cu  or  Zn)  have  been  ana¬ 
lyzed  using  neutron  diffraction  and  ab  initio  molec¬ 
ular  dynamics  (QMD)  modeling.  All  six  partial 
radial  distribution  functions  (PRDFs),  g/f  r),  have 
been  identified  for  each  alloy. 

2.  The  nearest-neighbor  mode  bond  lengths  are  short¬ 
er  than  those  in  competing  crystals.  A  noticeable 
shortening  of  Ca-Cu,  Mg-Cu,  and  Ca-Zn  bond  dis¬ 
tances  indicate  strong  interactions  between  these 
atom  pairs,  as  supported  by  the  calculated  elec¬ 
tronic  structure.  It  is  suggested  that  the  bond  short¬ 
ening  is  enabled  by  the  absence  of  long-range 
atomic  order,  which  lowers  the  free  energy  of  metal¬ 
lic  glasses  and  increases  GFA. 

3.  Pronounced  CSRO  near  TM  atoms,  chemical  short 
range  clustering  (CSRC)  near  Ca  atoms  and  a  neu¬ 
tral  environment  near  Mg  atoms  are  shown. 
Increasing  the  Cu  concentration  from  15  to  35 
atomic  percent  increases  the  total  coordination 
number  around  Ca,  CNCa,  from  13.6  to  15.0,  while 
CNMg  and  CNCu  remain  unchanged  at  « 12.5  and 
« 10.4,  respectively.  The  partial  coordination  num¬ 
bers  depend  on  alloy  composition,  so  that  the  num¬ 
ber  of  Ca  atoms  decreases,  the  number  of  Cu  atoms 
increases  and  the  number  of  Mg  atoms  is  almost 
constant  with  increasing  TM  concentration. 

4.  Voronoi  tessellation  show  that  many  types  of  coor¬ 
dination  polyhedra  are  present,  but  the  most  com¬ 
mon  are  Kasper  and  distorted  Kasper  polyhedra  at 
the  favorable  CNs.  The  fractions  and  distributions 
of  these  clusters  depend  on  alloy  composition.  Poly- 
tetrahedral-type  clusters  and  five-coordinated  verti¬ 
ces  dominate  in  the  amorphous  structures,  which 
indicates  that  tetrahedra  and  pentagonal  bi-pyramids 
are  the  main  building  blocks  in  these  amorphous 
alloys. 
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